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Abstract 

Positron lifetime spectra were measured in liquid water at temperatures between 0°C and 
50° C. The long lifetime of ort/io-positronium atoms (o-Ps) determined by electron pick-off in 
molecular substances decreases smoothly by 10% as the temperature is raised. This lifetime 
temperature dependence can be explained by combining the Ps-bubble model and the mixture 
state model of liquid water. 



I. INTRODUCTION 



Water has been studied extensively because it is a fundamental liquid with unique 
properties. For example, its boiling point under atmospheric pressure, 100°C, is unusually 
high for its molecular weight, 18. In addition, its volume decreases as it melts, unlike other 
liquids. Moreover, it has an anomalous density maximum at 4°C. These characteristics 
have been studied in detail using various methods. Consequently, many models have 
been suggested. Nevertheless, no definite mechanism has been established to explain 
these anomalous behaviors of water, even though each model can explain some specific 
phenomena. 

One sensitive method to investigate molecular states in solid and liquid forms is the 
positron annihilation lifetime (PALT) method. When positrons are injected into a liquid, 
some of them form positroniums with the surrounding electrons. The positronium lifetime 
is determined by the rate of annihilation between the positron in the positronium and 
surrounding electrons, also known as pick-off. As described later, according to the Ps- 
bubble model , the lifetime should be longer at higher temperatures in liquid water. 

However, Vertes et al. ^ showed that the positronium lifetime has a general decreasing 
trend for higher temperature, although it oscillates between 38°C to 60°C. 

We measured PALT in liquid water between 0°C and 50°C. This paper explains that the 
positron lifetime decreases smoothly as the temperature rises. This behavior is explainable 
by combining the Ps-bubble model and a two-state mixture model of water. 



II. POSITRON ANNIHILATION LIFETIME METHOD AND THE PS- 
BUBBLE MODEL 

This section explains PALT and the Ps-bubble model. 

Positrons injected in material can undergo three different processes, each with a dif- 
ferent lifetime. The shortest lifetime (ri) occurs by a formation and annihilation of para- 
positronium (p-Ps) in which the spins are anti-parallel. In that situation, the electron 
forming p-Ps with the positron is captured from the surrounding molecules. The lifetime 
of p-Ps in vacuum, ~ 0.1 ns, is too short to be affected by surrounding matter. 

The second lifetime {T2) is caused by positron annihilation without forming a bound 
state. 
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The third hfetime (rs) results from formation of an ort/io-positronium (o-Ps) with 
parallel spins. 

The o-Ps has an intrinsic lifetime of 142 ns in vacuum. However, in material, the 
lifetime is typically several nanoseconds because it annihilates with an electron in the 
surrounding molecules. As a result of this pick-off process, is sensitive to the electron- 
state of the surrounding substance. In liquid, o-Ps is generally considered to push out the 
surrounding molecules and form a Ps-bubble. The wave function of o-Ps exudes from the 
surface of Ps-bubble, which comprises electrons from surrounding substances. The overlap 
of o-Ps wave function and an electron wave function of the surroundings determines the 
pick-off rate. In the Ps-bubble model, the wave function of o-Ps is trapped inside a 
spherical infinite-depth well potential whose radius is the sum of the radius of the Ps- 
bubble and the exuding depth. Therefore, is a function of the Ps-bubble radius 

Nakanishi et al. introduced a semi-empirical correlation between r-? and the Ps-bubble 
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where is measured in nanoseconds, R is the Ps-bubble radius, and AR is the exuding 
depth of the o-Ps wave function into the surrounding electron wave functions. In this 
model, R is determined by the balance between the zero-point energy of o-Ps and the 
surface tension of the surrounding substance. The balance is represented as 

-^(E + AnR'^)=0, (2) 

where E is the zero point energy of o-Ps and 7 is the surface tension of bulk water. E is 
expressed as 

E ~ = (3) 

Arue Ame{R + ARy' ^ ' 

where me is the electron mass and k is the momentum of o-Ps. 

For the case of water, 7 decreases as the temperature rises; thereby the bubble becomes 
larger. By this argument alone, is longer for higher temperatures. 
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FIG. 1: A schematic of the apparatus, a. BaF2 scintillators with photomultiplier tubes, b. 
Sample vial, c. Satellite water bath (60 mm x 60 mm x 100 mm (height)), d. Main water bath 
(300 mm x 500 mm x 160 mm (height)), e. Needle for pressure equilibrium, f. K-thermocouple, 
g. Positron source. 

III. EXPERIMENT 

A. Measurement of positron lifetime 

We employed the positron annihilation lifetime technique using ^^Na as a positron 
source. ^^Na is put in the sample center. Positrons emitted into the sample are annihilated 
through the process mentioned in the previous section. When ^^Na emits a positron, it 
simultaneously emits a 1275 keV photon. When a positron is annihilated with an electron, 
two 511 keV photons are emitted. In light of those facts, we measured the time difference 
between the emission of 1275 keV photon and the emission of two 511 keV photons. 

We used the apparatus shown in Fig. [T] Powder of 1.6 MBq ^^NaCl wrapped between 
8 /im thick polyimide films was used as a positron source. The polyimide film is known 
to have no rs component The source was held inside a 20 mm diameter 25 ml glass 
vial. The vial was then filled with de-gassed distilled water. We estimate that 4.4% of the 
positrons were absorbed in polyimide film; more than 95.5% of the positrons are absorbed 



within 0.5 mm of water ^. The vial was submerged in a water bath to maintain the sample 
water temperature within ±0.1°C. The sample water pressure was maintained at 1 atm 
by inserting a 0.9 mm inner diameter needle into the vial. 

The 1275 keV photon (start signal) was detected using a barium fluoride (BaF2) scin- 
tillator (38.1 mm diameter, 25.4 mm thickness) with a photomultiplier tube (H3378-51; 
Hamamatsu Photonics, K.K.) and one 511 keV photon (stop signal) was detected by an- 
other BaF2 detector. The two BaF2 detectors viewed the vial at a 90° opening angle; they 
were held 27 mm away from the center of the sample water. 

The time difference between the start and stop signals was measured with a time-to- 
amplitude convertor and a multi-channel analyzer. The time resolution of this system is 
~ 280 ps. We measured the positron lifetime at 14 temperatures between 0°C and 50°C. 
In most cases, we performed four runs at each temperature point. We used a new water 
sample for each run. Each run collected one million events in 1500 s. 



B. Analysis 

A typical lifetime spectrum acquired in our experiment is shown in Fig. |21 
To extract lifetimes, the spectrum was fltted for a model function, 

3 3 

S(t) = J2- [ e-''/^^git-f)9it' -to)dt' + kJ2- [ e'' /^'g{t - t')9ito - t')dt' + C, (4) 
i=i ^ i=i ^ 

where i is the index for different lifetimes, Jj is intensity, Tj is the lifetime, and g{t — t') 
is the Gaussian describing the time resolution. Also, 6{t) represents a step function in 
which 6{t) = 1 (t > 0) and (t < 0), to is the delayed timing of the start signal, and C is 
a constant background rate. The second term is added to accommodate events in which 
a 511 keV photon gave the start signal and a 1275 keV photon gave the stop signal. Ten 
)arameters, /i,2,3) ''"1,2,3) ^05 k, C, and the sigma of the Gaussian were fltted using MINUIT 
^ fltting code, which is commonly used in analyses of high-energy physics. A fltting 
result achieved using MINUIT is shown as a solid line in Fig. El 

We used POSITRONFIT jl^, which is a popular code in positron lifetime analysis, 
only to crosscheck our results because POSITRONFIT requires the time resolution as an 



^ Using the linear absorption coefficient (a), a = Ad / E^g^-^{cm where d is the density of absorber 
and -©max the maximum energy of a positron in MeV. 
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FIG. 2: A typical positron lifetime spectrum in liquid water (bottom) and the data-to-fit ratios 
(top). Data are shown in dots. A fitting result by MINUIT is shown in a solid line in the 
spectrum. 

input parameter. The results are sensitive to the lower boundary of the fitting range. 
Both codes gave mutually consistent results. The mean value of the difference between 
both codes taken over all measurements was 0.008 ns; their standard deviation was 0.016 
ns. 



IV. RESULTS 

Temperature dependence of is shown in Fig. El The decreases smoothly by 10% 
as the temperature is raised from 0°C to 50°C. Our ts at 20°C, 1.839 ± 0.015 ns, agrees 
with the measurement of Mogensen, 1.85 ns |ll(. 

Behaviors of other lifetime parameters are shown in Fig. 0] We confirmed that the 
behavior of is not caused by a change in other lifetime parameters. To do so, we 
produced artificial spectra that represent a sum of three exponentials, Yl^=i 6xp(t/rj) / Tj, 
convoluted by a Gaussian time resolution. The number of events in each 0.0495 ns time 
bin was made to follow a Poisson distribution. We produced such artificial data samples 
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FIG. 3: T3 in water as a function of temperature. Dots show data. Vertical lines show errors. 
The solid line shows a fitting result by Ps-bubble model combined with the two-state model. 
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FIG. 4: Temperature dependencies of lifetimes, ri and T2 in water and intensities of three 
lifetime components. Dots show data and vertical lines indicate the fitting error. Solid lines are 
results of linear regression, a, b: lifetime component ti,T2, respectively, c: intensities of three 
components. 



with different Xj and Jj, and fitted for them. Figure El shows the fitted r-^ as a function 
of the varied parameters. Solid hues show the linear fit to the points. First, the fitted 
Ts agrees with the input within 0.01 ± 0.01 ns, as shown in Fig. El^c). Second, the 
fitted r3 does not depend on ri, T2, I2, and I3. Horizontal bars show the variations and 



7 



errors of parameters in the measured temperature range. The maximum deviation of 
Ts is calculated by multiplying the slope of the linear fit and the quadratic sum of the 
variation and the error of the other parameters. Deviations resulting from ri, T2, I2, and 
I3 are 0.000 ± 0.004 ns, 0.0014 ± 0.0008 ns, 0.0015 ± 0.0029 ns, and 0.001 ± 0.003 ns, 
respectively. The only apparent dependence is on Ji, as shown in Fig. El^d); its effect 
on Ts is 0.009 ± 0.002 ns. The total effect on T3 that is attributable to the error and 
swing of other parameters in the temperature range is 0.011 ns. Consequently, we can 
conclude that the temperature dependence of is not caused by changes in other lifetime 
parameters. 




12 14 16 18 20 22 24 

h (%) 



20 21 22 23 24 25 
/3 (%) 



FIG. 5: Fitted as a function of other parameters studied with artificial data samples: (a) ri, 
(b) T2, (c) r3, (d) /i, (e) I2, and (f) I3. Horizontal solid bars show averages of the standard 
deviations of parameters at each temperature point; dashed bars show the maximum changes of 
parameters between and 50°C. 



V. DISCUSSION 

The behavior of our cannot be explained by the original Ps-bubble model. For that 
reason, we apply the two-state mixture model to the Ps-bubble model. Vedamuthu et al. 
showed that the two-state mixture model can reproduce density of H2O in the temper- 
ature range between -30 and +70°C with very high precision [l^]. In this model, liquid 
water comprises two dynamically inter-converting mixed micro-domains whose bonding 
characteristics are similar to those of ice \h (lower density) and ice II (higher density) 
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. The total specific volume of liquid water is given as 

V{T, P) = [l- /ii(T, P)]V^{T, P) + /ii(T, P)V^i{T, P), (5) 



where I (II) indicates lower (higher) density bonding, /jj is the mass fraction of the higher 
density bonding type, and ^(^n) is the specific volume of the lower (higher) density type 
bond. 

We modified the Ps-bubble model by introducing the two-state mixture model with 
the following assumptions. 

1. Positronium forms a bubble (Ps-bubble) with a radius R in water as the original 
Ps-bubble model. 

2. Water consists of two different molecular bond types, I and II; the ratio /jj is a 
function of the temperature as given by Vedamuthu et al. |12]. 

3. The pick-off rate is determined by an overlap between the o-Ps wave function and 
the electron wave function of the water; it can be parameterized by exuding depths, 
A_Rj and A_Rjj for each bonding type. 

4. The bubble radius, i?, is given by a macroscopic surface tension 7, which is deter- 
mined by the temperature jl^, Ai?j and Ai?jj . 

The modified equation (jT)) is: 

R 1 / 2tiR n ^ n -1 



rr i?+(i-/n)Ai?i + /nAi?n '^27r''''(i? + (i-/n)Ai?i + /nAi?n)jJ 

Figure El shows the fitted using Eq. © as a function of temperature. The fitting 
parameters are Ai?j and Ai?jj. The reduced is 1.7. The modified Ps-bubble model 
represent the temperature dependence of rs well. 

The fitted exuding depths are Ai?j = 0.130 ± 0.005 nm and Ai?jj = 0.218 ± 0.003 
nm. In the Ps-bubble model, Ai? should correlate with the van der Waals radius of the 
surrounding atoms. In fact, Ai?, which is determined by fitting data of well-characterized 
small-pore materials such as zeolites for Eq. (Q), is 0.166 nm 0, 15|, which agrees with 
the 0.166 nm average of van der Waals radii of atoms of Si04, the main component of 
zeolite. In our result, Ai?j = 0.130 nm is close to the van der Waals radius of oxygen 
(0.155 nm) and hydrogen (0.120 nm). Contrarily, Ai?jj is too large for the van der Waals 
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FIG. 6: Temperature dependence of the Ps-bubble radius, as calculated with 7 and fitted Ai?j 
and Ai?jj. 



radius. This indicates that the higher density state has a larger contribution to the pick- 
off process than their own depth of the electron wave function. This constitutes evidence 
that the higher density state, which has a bending hydrogen bond is more active 
than the lower density state. 

Next, the fitted R supports the Ps-bubble idea. As shown in Fig. El the radius of the 
Ps-bubble, R, based on our fit, is about 0.3 nm. Liquid water has intrinsic vacancies in 
its structure while its structure is changing rapidly. Those vacancies produce a hexagonal 
structure having an oxygen atom at each vertex connected by hydrogen bonds. Consider- 
ing the distance between the nearest O- ■ -O, 0.28 nm Ig], and the van der Waals radius 
of water, the vacancy is smaller than a Ps-bubble. Therefore o-Ps should spread water 
molecules apart to exist in water. This fact supports the need for a balance between the 
zero-point of o-Ps energy and the surface tension of water. 

In addition, our model shows that R has a minimum at 8°C. The radius of Ps-bubble 
at 50°C is 1.009 times the radius at 8°C. This size relationship is consistent with the fact 
that the distance between the nearest neighbors of O- ■ -O pair at 50°C is 1.011 times the 
distance at 4°C Q. 
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VI. CONCLUSION 



We precisely measured the temperature dependence of the long lifetime of a positron 
(ts) in water at temperatures between and 50°C. The T3 decreases smoothly as temper- 
ature rises; it is 1.839 ±0.015 ns at 20°C. This fact is explained by combining two models 
in which the water consists of two types of molecular bonds, and in which positroniums 
push those bonds apart to form Ps-bubbles. We also found that T3 is sensitive to the 
electron states of the two types of molecular bonds. 
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